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TRACE ESTIMATES FOR UNIMODAL LEVY PROCESSES 


K. BOGDAN AND B. A. SIUDEJA 


Abstract. We give two-term small-time approximation for the trace of the Dirich- 
let heat kernel of bounded smooth domain for unimodal Levy processes satisfying 
the weak scaling conditions. 


1. Introduction 

A two-term small-time uniform approximation for the trace of the transition density 
of the Wiener process killed off bounded R-smooth domain D C i.e. the classical 
Dirichlet heat kernel, was obtained by van den Berg [16]. The first term of the 
approximation is proportional to the domain’s volume \D\ and the second-to the 
surface measure \dD\ of the boundary, with explicit coefficient depending on time. 
Asymptotic non-uniform expansions of the trace of the heat kernel were given earlier 
in [11], see the discussion in [16]. 

Banuelos and Kulczycki [1] obtained a uniform two-term approximation for the 
isotropic a-stable Levy processes. The closely related case of the relativistic a-stable 
Levy processes was resolved by Banuelos, Mijena and Nane [3]. A similar two-term 
approximation for Lipschitz domains was given for the Wiener process by Brown [8], 
and for the isotropic o-stable Levy processes-by Banuelos, Kulczycki and Siudeja [2], 
Park and Song [12] obtained a two-term small-time approximation of the trace for 
the relativistic a-stable Levy processes on Lipschitz domains, and gave an explicit 
power expansion of the first term. 

In this work we investigate those Levy processes X t in M. d , where d > 2, which are 
unimodal and satisfy the so-called weak lower and upper scaling conditions, denoted 
WLSC and WUSC respectively, of orders strictly between 0 and 2 (see Section 2 for 
details). The isotropic stable and relativistic Levy processes are included as special 
cases but at present the orders of the lower and upper scalings may differ. For bounded 
R-smooth open sets D C (also called C 1,1 open sets in the literature) our main 
result gives a two-term small-time approximation of the trace of the corresponding 
Dirichlet heat kernel. For instance we resolve sums of independent isotropic stable 
Levy processes with different indexes. 
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In what follows we let -0 be the Levy-Khintchine exponent and pt(x) be the tran¬ 
sition density of X t . We consider 

td — {t > 0 : X t 0 D}, 

the hrst time that X t exits D. For t > 0 and x, y G we dehne the heat remainder 
r D (t,x,y) = E x [r D < t,p t _ TD (X(r D ) - y )\. (1) 

The Dirichlet heat kernel for X t is given by the Hunt formula: 

PD(t,x,y) = p t (y - x) - r D (t,x,y), (2) 

and the trace of X t on D is 

tr (t, D) = j po(t,x,x)dx, t > 0. (3) 

We denote HI = {(or,..., Xd) G : X\ > 0}, a half-space, and for t > 0 we let 

POO 

C m (t) = / r m (t, (q, 0, • • • , 0), (q, 0, • • ■ , 0 ))dq. 

Jo 

For instance, C'n(t) = ct~ d ^ a+1 ^ a for the isotropic a-stable Levy process [1]. Here is 
our main result (a stronger statement is given as Theorem 3.1 in Section 3). 

Theorem 1.1. If bounded open set D C R d is R-smooth, WLSC and WUSC hold 
for if, and t —> 0, then tr(t, D) equals p t (0)\D\ — Cm(t)\dD\ plus lower order terms. 

Heuristically, if x G D and t > 0 is small, then r D (t, x, x) is small and so po(t, x, x) 
is close to p R d(t, x,x) =p t (0). Therefore the first approximation to tr(t, D) isp t (0)|D|. 
The second term in Theorem 1.1, Cm(t)\dD\, approximates J D r D (t,x,x)dx. As we 
shall see, r£>(t,x,x ) depends primarily on the distance of x from dD. It is here 
that the A-smoothness of D plays a role by allowing for an asymptotic coefficient 
independent of D, that is C'n(t). In view of the definition of C^t), the appearance 
of \dD\ in the second term of the approximation of the trace is natural. 

In some cases, including the relativistic stable Levy process, explicit expansions of 
Pt(0) can be given [12, Lemma 3.2], In more general situations p t (0), Ce(t) and the 
bounds for the error terms cannot be entirely explicit but Lemma 2.7 and Theorem 3.1 
below provide a satisfactory formulation. 

Technically we only need to estimate f D rr>(t, x, x)dx to prove Theorem 1.1. In this 
connection we note that sharp global estimates for p£>(t,x,y ) were recently obtained 
by Bogdan, Grzywny and Ryznar [6], but these estimates do not easily translate into 
sharp estimates of ro(t,x,y). Namely, if po(t,x,y) is only known to be proportional 

Pt(y ~ x ), then essential further work is needed to accurately estimate r D (t,x,y). 

The paper is composed as follows. In Section 2 we give preliminaries on unimodal 
Levy processes with scaling, their heat kernel, Green function and Poisson kernel 
for .R-smooth open sets. In Section 3 we prove Theorem 3.1, a stronger and more 
detailed variant of Theorem 1.1. The most technical step of the proof of Theorem 3.1 
is given separately in Section 4. 
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We remark in passing that the trace can also be studied and interpreted within the 
spectral theory of the corresponding semigroup given by the integral kernel pr> [1], 
In view toward further research we note that sharp pointwise estimates of r£)(t,x,y) 
complementing [6] would be of considerable interest. We also note that two-term 
approximations of the trace of the heat kernel of general unimodal Levy processes are 
open for Lipschitz domains. 

Acknowledgments. We thank Tornasz Grzywny for very helpful discussions and 
suggestions on the manuscript. 


2. Preliminaries 

2.1. Unimodality. A Borel measure on R d is called isotropic unimodal, in short: 
unimodal, if on R d \ {0} it is absolutely continuous with respect to the Lebesgue 
measure and has a radially nonincreasing, in particular rotationally invariant, or 
isotropic density function. Recall that Levy measure is an arbitrary Borel measure 
concentrated on R d \ {0} and such that 

/ (|x| 2 A l) u(dx) < oo. 

J R d 

In what follows we assume that v is a unimodal Levy measure and define 

^(0 = [ (1-cos (£,x))v(dx), £eR d , (4) 

jR d 

the Levy-Khintchine exponent. It is a radial function, and we often let 0(r) = 0(£), 
where £ G R c/ and r = |£| > 0. The same convention applies to all radial functions. 
The (radially nonincreasing) density function of the unimodal Levy measure v will 
also be denoted by u, so u(dx) = u(x)dx and u(x) = v(\x\). We point out that for 
A > 1 and r > 0, ijj(Xr) > Tr~ 2/ ip(r) and 0(Ar) < 7r _2 A 2, 0(r) [5, Section 4], More 
restrictive inequalities of this type define what are called the weak scaling conditions, 
see Section 2.2. 

We consider the pure-jump Levy process A" = (X t , t > 0) on R d [13], in short: X t , 
determined by the Levy-Khintchine formula 

Ee i<£,x t ) =e -m) = f e i{ S’ x) p t (dx). 

Jm d 

The process is (isotropic) unimodal, meaning that all its one-dimensional distributions 
Pt(dx) are (isotropic) unimodal; in fact the unimodality of v is also necessary for the 
unimodality of X t [17]. In what follows we always assume that ^ is unbounded, 
equivalently that u(R d ) = oo. In other words X t below is not a compound Poisson 
process. Clearly, -0(0) = 0 and ^{u) > 0 for u > 0. By [6, Lemma 1.1], pt(dx) 
have bounded, in fact smooth density functions Pt(x) for all t > 0 if and only if the 
following Hartman-Wintner condition holds, 


lim 0(£)/ln|£| = oo. 
ki-*» 


( 5 ) 
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Let V be the renewal function of the corresponding ladder-height process of the 
first coordinate of X t . Namely we consider X^\ the first coordinate process of X t , 
its running maximum M t := sup 0<s<i Xs 1 ^ and the local time L t of M t — x[ 1 ^ at 0 so 
normalized that its inverse function Z0 1 is a standard 1/2-stable subordinator. The 
resulting ladder-height process rj(t) := X^’^Lj 1 ) is a subordinator with the Laplace 
exponent 


k{u) = -logEe-^ (1) = 


exp 


log 0(0 

1 + C 2 


dC 


u > 0, 


and V(x) is defined as the accumulated potential of r/: 


V(x) = E / 1[ 0 )X ]{r]t)dt, x>0. 

Jo 

For x < 0 we let V(x) = 0. For instance, if 0(£) = |£| a with a G (0,2), then 
V(x) = x" /2 [15, Example 3.7]. Silverstein studied V and V' as g and 0 in [14, (1.8) 
and Theorem 2], The Laplace transform of V is 

1 

/ V(x)e 


c dx = 


u > 0. 


UK(U 


We 


( 6 ) 


The function V is continuous and strictly increasing from [0, oo) onto [0, oo 
have lim^oo V(r) = oo. Also, V is subadditive: 

V(x + y) < V(x) + V(y), x, y G M. 

For a more detailed discussion of V we refer the reader to [4] and [14]. 

In estimates we can use V and 0 interchangeably because by [6, Lemma 1.2], 

V{r) pa [0(l/r)] -1 / 2 , r > 0. (7) 

The above means that there is a constant, i.e. a number C G (0, oo), such that for all 
r > 0 we have C _1 E(r) < [0(l/r)] _1 ^ 2 < CV (r). In fact in (7) we have C = C(d), 
meaning that C may be so chosen to depend only on the dimension, see ibid. Similar 
notational conventions are used throughout the paper. To give full justice to V, the 
function is absolutely crucial in the proofs of [4], a paper leading to [6]. By (6), 


-eV(r) < V(er) < V(r), 


0 < £ < 1, 0 < r < oo. 


( 8 ) 


2.2. Scaling. We shall assume relative power-type behaviors of 0(r) at infinity. 
Namely we say that 0 satisfies the weak lower scaling condition at infinity (WLSC) 
if there are numbers a > 0, 6 G [0, oo) and C_ G (0,1], such that 

0(A r) > C\ -0(r) for A > 1, r > 6. 

Put differently and more explicitly, 0(r)/r- is almost increasing on (0, oo), i.e. 


VKf) > c Mj) 


rr* U 


if 


s > r > 6. 


In short we write 0 G WLSC (a, 9, C ), 0 G WLSC(a, 0), 0 G WLSC(a) or 0 G WLSC, 
depending on how specific we wish to be about the constants. If 0 GWLSC(a, 0), then 
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we say that 0 satisfies the global weak lower scaling condition (global WLSC) if 0 = 0. 
If 0 > 0, then we can emphasize this by calling the scaling local at infinity. We always 
assume that 0^0, therefore in view of -0 gWLSC we have the Hartman-Wintner 
condition (5) satisfied, and so B0 3 x H > pt(x ) is smooth for each t > 0. 

Similarly, the weak upper scaling condition at infinity (WUSC) means that there 
are numbers a < 2, 0 > 0 and C£ [1, oo) such that 

0(Ar) < CA“0(r) for A > 1, r > 0. 

In short, 0 £ WUSC 000,(0) or 0 £ WUSC. Global WUSC is WUSC(a,0), etc. 

We call a, 0, C_, a, 9 , C the scaling characteristics of 0. As pointed out in [6, 
Remark 1.4], by inflating C and C we can replace 9 with 0/2 and 6 by 0/2 in the scal¬ 
ings, therefore we can always choose the same, arbitrarily small value 0 = 0 = 0 > 0 
in both local scalings WLSC and WUSC, if they hold at all. The scalings character¬ 
ize the so-called common bounds for p t {x) [5, Theorem 21 and Theorem 26], and so 
they are natural conditions on 0 in the unimodal setting. The reader may also find 
in [5] many examples of Levy-Khintchine exponents which satisfy WLSC or WUSC. 
For instance 0(£) = |£|“, the Levy-Khintchine exponent of the isotropic a-stable 
Levy process in 10 with a £ (0, 2), satisfies WLSC(a, 0,1) and WUSC(a, 0,1). The 
characteristic exponent 0(£) = (1 + \f\ 2 ) a ^ 2 — 1 of the relativistic ct-stable Levy pro¬ 
cess with a £ (0,2) satisfies WLSC(cpO) and WUSC(a, 1). Other examples include 
0(0 = |£| ai + |<el“ 2 e WLSC(«i, 0, l)OWUSC(a 2 , 0,1), where 0 < or < a 2 < 2, etc. If 
0(r) is a-regularly varying at infinity and 0 < a < 2, then 0 £WLSC(a)flWUSC(a), 
with any 0<a<a<a<2. The connection of the scalings to the so-called 
Matuszewska indices of 0(r) is explained in [5, Remark 2 and Section 4], 

If 0 £WLSC(a, 0), then by (7) (or see [6, (1.8)]) we get the following scaling at 0: 

V{er ) < Ce~ /2 V (r), 0 < e < 1, 0<r< 1/0. (9) 

Here the range is 0 < r < oo if the lower scaling of 0 is global, in agreement with (9) 

and the convention 1/0 = oo. If 0 £WUSC(a, 0), then, similarly, 

V{er) > C£“ /2 U(r), 0 < e < 1, 0 < r < 1/0. (10) 

We shall need U -1 , the inverse function of V on [0, oo). We let 

T(t) = V-\Vi), t > 0. (11) 

Put differently, [V(T(t))] 2 = t. For instance, T(t) = t l ^ a for the isotropic a-stable 
Levy process. The functions V and T allow us to handle intrinsic difficulties which 
hampered extensions of [16, 1, 3, 12] to general unimodal Levy processes, namely the 
lack of explicit formulas and estimates for the involved potential-theoretic objects. 

We note that T(t) < a if and only if t < V 2 (a), wherever a, t > 0. The scaling 
properties of T at zero reflect those of 0 (at infinity) as follows. 

Lemma 2.1. If (9) holds, 0 < £ < 1 and 0 < t < U(l/0) 2 , then T{et ) > ce 1// -T(t). 
If (10) holds, 0 < £ < 1 and 0 <t< U(l/0) 2 , then T(et ) < c£ l ^ a T{t). 
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Proof. To prove the first assertion we note that T is increasing. If 0 < t < V{1 /9) 2 , 
and 0 < e < 1, then T(t ) < 1/9 and T(et)/T(t) < 1. By (9), 

a/2 




(T(st)Y 

V T(t) ) 


□ 


V(T(t)) 

as needed. The proof of the second inequality is analogous but uses (10). 

By (8) and the proof of Lemma 2.1 we always have 

T(et) < cy/eTit), 0 < e < 1, 0<r<oo. (12) 

In what follows we always assume that v is an infinite unimodal Levy measure on 
W l with d > 2 and the Levy-Khintchine exponent defined by (4) satisfies 

if G WLSC(a, 9) n WUSC(a, 9 ), 

where 0 < a < a < 2, and 9 > 0. Many partial results below need less assumptions 
but for simplicity of presentation we leave such observations to the interested reader. 

Definition. We say that (H) holds if for every r > 0 there is H r > 1 such that 

V(z) — V ( y ) < H r V\x)(z — y) whenever 0 < x < y < z < 5x < 5r. 

We say that (H*) holds if Hoo := sup r>0 ihr < oo. 

We may and do chose H r nondecreasing in r. By [4, Section 7.1], (H) always holds 
in our setting because if satisfies WLSC and WUSC. If GWLSC(a, 0)flWUSC(a, 0), 
then (H*) even holds. 


2.3. Heat kernel. By [6, Lemma 1.3], there is a C\ — C\(d ) such that 

Pt{x) - Cl WVY\^y xeR d \{0}, (13) 

hence [5, (15)], 

V{X) S ° l V 2 (\x\)\x\ d ' X * °' (14) 

Since if eWLSC(a,6 1 ), by [6, Lemma 1.5] we have 

p t (x) < cT~ d (t), t<V 2 (9~ l ), xeR d . (15) 

We now discuss the heat remainder and the heat kernel of open sets D C W l . As 
usual, 0 < r D (t,x,y ) < p t (x — y ). Indeed, one directly checks that [0 ,t) 3 s H>• 
Y s = p(t — s,X s ,y) is a P x -mart ingale for each x,y G M d . The martingale almost 
surely converges to 0 as s —* t, and we let Y t — 0. By optional stopping, quasi¬ 
left continuity of X and Fatou’s lemma, for every stopping time T < t we have 
E x Y t < E X Y 0 = p(t,x,y). The inequality r D (t,x,y ) < p t (x — y) follows by taking 
T = r D A t. The next result is a consequence of the strong Markov property of X t . 

Lemma 2.2. Consider open sets D C F C M rf . For all t > 0 and x, y G M. d , 

p F (t, x, y) - p D (t,x,y ) = E y [r D < t,X(r D ) G F\D; p F (t - t d ,X(t d ),x)] . 






TRACE ESTIMATES FOR LEVY PROCESSES 


7 


Proof. We repeat verbatim the proof of [1, Proposition 2.3]. □ 

Here is a well-known Ikeda-Watanabe formula for the joint distribution of X{td ) 
and td, see [10, Proposition 2.5] or [7, (27)] for proof. 

Lemma 2.3. Let D C M. d be open. For x G D, t 2 > t\ > 0 and A C ( D) c , 

Pd(s, x, y)ds / u(y — z)dzdy. 

J A 

We denote 8d(x) := dist(x, D c ), x G R d . 

Lemma 2.4. We have 


¥ x (X(td) G A, ti < To < t 2 ) — I f 

J D J t\ 


and 


r D (t,x,y ) < C i 


r D (t,x,y) < CT(t) 
t 


x,y e 


(16) 


(17) 


v 2 (M*))W 

Proof. Since 0 G WLSC(a, 6 1 , C), we have (15), which yields (16). By (1), (13), and 
symmetry, 

ro{t,x,y) = r D (t,y,x) < E y t d < t; Ci- 


t — t d 


V*{\X{t d ) ~ x\)\X{t d ) - x\\ 

Since |X(td) — x\ < 5 d(x) and V is increasing, we obtain (17). □ 

Recall that El is a half-space and Ce(£) is defined immediately before Theorem 1.1. 
Lemma 2.5. IfT(t) < 1/9, then Cjg(f) < cT(t)~ d+1 . 

Proof. Denote r(t, q) = rm(t, (q, 0, • • • , 0), (q, 0, ■ ■ ■ , 0)). By (17) and (9), 


r(t, q)dq < c 


"OO t/2 


JT(t) 

Using (16) we get 


'm 


V\T(t)) 

V 2 (q)q d 


dq < c 


r»oo 

l , 00^= 

' T(t) q - 


1 -d 


rT(t) 


rT(t) 


r(t,q)dq<c I T(t ) d dq = cT{t ) 


\l-d 


□ 


To obtain a lower bound for Cn(t) we shall use the existing heat kernel estimates 
for geometrically regular domains. Recall that open set D C satisfies the inner 
(outer) ball condition at scale R > 0 if for every Q G dD there is a ball B(x', R) C D 
(a ball B(x",R) C D c ) such that Q G dB{x',R ) (Q G dB(x",R), respectively). An 
open set D is A-smooth if it satisfies both the inner and the outer ball conditions at 
some scale R > 0. We call B(x',R) and B(x”,R) the inner ball and the outer ball, 
respectively. 

In the next lemma we collect a number of results from [6]. For brevity in what 
follows we sometimes write T = T(t), where t > 0 is given. 
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Lemma 2.6. Let open D C W l satisfy the outer ball condition at scale R < 1/9. 
There is a constant c such that for TV \x — y\ < 1/6, 


Pd (t, x, y) < c 


( V{Sp{x)) 
\V(T A R) 



( V(S D (y)) 

\V(TAR) 



(r~ d A 


V 2 (T) 
x — y\ d V 2 (\x 



Proof. We have (H). We note that \ft, — V ( T ) and use the second part of [6, Corollary 
2.4], We need to justify that the quotient Hr/J 4 (R) is bounded, where Hr is the 
constant from (H) and J(R) = info< r <_R v(B(0, r) c )V 2 (r). To this end we observe 
that Hr is increasing, and J(R) is nonincreasing, hence we get an upper bound for 
this quotient by replacing R with 1/9. If 9 = 0, which we also allow, then by [4, 
Proposition 5.2, Lemma 7.2 and 7.3] the quotient is bounded as a function of R. By 
[6, Lemma 1.6] with r = 1/2, we also have p t / 2 (0) < cT~ d {t). □ 

Lemma 2.7. We have Cm{t) ~ T(t)~ d+1 ~ p t (0)T(t) as t —> 0. 


Proof. By Lemma 2.6 and (2) there is s > 0 such that r(t, q) > ^Pt(O) if V(q) < £\/i. 
Since -0 GWUSC, by scaling of V there is c > 0 such that for 0 < q < cT(t) the 
condition is satisfied and we have 

fcT(t) i 

J o r(t, q)dq > - J T(t)~ d dq = -T{tf~ d . 

By WUSC and WLSC we have p t ( 0) « T[f)~ d , see [5, (23)]. □ 


2.4. Green function. For M > 0, the truncated Green function of D is defined as 


-M 


G D (x,y)= p D {t,x,y)dt , x,yeR. 

Jo 

The Green function of D is 

DOG 

Gp(x,y)= p D {t, x, y)dt — G™(x, y). 

Jo 

Lemma 2.8. Let open D C M d satisfy the outer ball condition at scale R < 1/9 
x, y G and \x — y\ < 1/9. Let M = V 2 (R). Then 

V{6 D {y))V(S D (x)) 


G d (x,y) < c 


and 


G d (x,y) < c 


\x — y\ d 

V{5 D (y))V(\x-y\) 
\x — y\ d 


(18) 


(19) 


Furthermore, if d > 2 or WUSC(a,0) holds, then (18) and (19) even hold for M = 
V 2 {l/9), including the case of global WLSC (M = oo). 

Proof. Assuming T < R A \x — y\, by Lemma 2.6 we get 

V{T A 5 D {x)) 


V 2 {\x-y\)\x-y\ d ' 


Pp{t,x,y ) < cV(S D (y)) 
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hence 


-V 2 0 - 2 / 1 AR) V(S D (x)) 

p D (t, x, y)dt < c 


"V 2 (\x—y\AR) 


V(T A S D (x))dt 


< c 

< c 


V 2 (\x - y\)\x - y\ d J 0 
V{5 d {x))V\\x -y\t\ R)V{\x - y\ A M^)) 
\x — y\ d V 2 (\x — z\) 
V(5 d (x))V(\x — 2/| A S d (x)) 


\x — y\ d 

This establishes (19) and (18) for small times. Then, 


-V 2 (i?) 


•> V 2 (\x-y\) 

By WUSC and Lemma 2.1, 


With this in mind we obtain 

rV\R) T -d^ 


p D (t,x,y)dt < cV(5 d (x)) 


,V 2 (R) rp-d 


' V 2 (\x-y\) 


T~ d (t) 

7T 


dt. 


1 ce 1 /" 


T(t) ~ T(et) 


'V 2 (\x-y\) 


Vt 


dt < c 


= c 


V 2d/a (\x-y\) 


-dt 


JvHw-,n *‘ VB+1/2 T d (v\\x - j/D) 

V^(\x-y\) [v2{]x _ 


\x — y\ d 

where the integral converges, because d/a + 1/2 > 1 (recall that a < 2). We thus get 
(19). To finish the proof of (18) we note that 


»v 2 (r) 


JV 2 {\x—y\) 

and we proceed as before. 


p D {t,x,y)dt < cV(8 D (x))V(8 D (y )) 


rV 2 (R) T -d( t j 
V 2 (\x—y\) t 


dt, 


□ 


2.5. Poisson kernel. For M > 0, the truncated Poisson kernel is defined as 
Kd(x, z) = f Gp (x, y)u(y - z)dy , x G D, z e D c . 

J D 

Lemma 2.9. Let open D C R d satisfy the outer ball condition at scale R. // diam(ZXJ 
{z}) < 1/6, then 


k v(r 2 )/2, x < c 

D ’ — V(5d{z)) \x — z^' 


x e D, z e D c 


Proof. The previous lemma gives an estimate for G X D ( ' R \ and the Levy measure is 
controlled by (14). Thus, 

V{\x-y\)/\V(5 D (y)) 


K 


V(R 2 )/2 

D { 


x,z) < cV(8 d (x)) 


Id \x - y\ d \y - z\ d V 2 (\y - z\) 


dy. 
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Note that \x — y\ > \x — z\/2 or \y — z\ > \x — z\/2. Furthermore, if \x 
then \x — y\ > \x — z\/2. Therefore, it is enough to verify that 


/ := 
II : = 


and 


l D \y ~ z\ d V 2 (\y - z\) V(5 d (z )) 
/' V{\x-y\) 


I Dn{\x-y\<\y-z\} \x ~ y\ d V 2 (\y ~ Z |) 

Considering / we note that Sn(y) < \y — z\, hence 

r \ y _ z \-d 


dy < 


C 


I < 


'\y — z\>Sjy(z) ^ 1 ) 


dy < c 


V{Sd{z)Y 
l ' d dr 


'Sd(z) 


rV (r) 


Using the scaling (9) we get 


y I > \y-z\, 


J < C r°° 1 5 D (z) \- /2 dr c 

- V(5d(z)) Js d ( z )\ r J r V(S D (z))' 

To verify the estimate for II we also use the scaling properties of V. For y £ D we 
have \y — z\ < 1/9 , hence 


II < c 


< 


'\x-y\<\y-z\ 

C [ 


\x - y\ 
\y- A 

\y~ z \ 


a/2 


dy 


V(6 d {z)) Jo 


x-y\ d V(\y - z |) 

c 2 


r \ d r 


\y-z\J r V(5 D (z))a' 


□ 


3. Proof of the main result 


For the convenience of the reader in the following statement we repeat our standing 
assumptions; see also the definition of V in Section 2.1 and that of T in (11). 

Theorem 3.1. Let v be an infinite unimodal Levy measure on with d > 2, and 
let the Levy-Khintchine exponent (4) satisfy G WLSC(o, 6) D WUSC(a, 6), where 
0 < a < a < 2 and 6 > 0. Let open bounded set D C be R-smooth with 0 < R < 
1/9. There is a constant cy depending only on v and 9 such that if 0 < t < U 2 (6 )_1 ), 
or T(t) < 1/9, then the trace (3) of the Dirichlet heat kernel (2) satisfies 

T(t ) 2 


tr {t,D) - \D\p t (0) + \dD\C m (t) < ce\D\p t (0) 


R 2 


( 20 ) 


If 9 = 0, then (20) holds for all t > 0. 


Recall that Lemma 2.7 asserts that Ch(£) ~ pt(0)T(t) and p t ( 0) ~ T(t) d as t —> 0, 
so the approximation of the trace in Theorem 3.1 is given in terms of powers of T(t). 


Proof of Theorem 1.1. The result is a direct consequence of (15), Lemma 2.7 and 
Theorem 3.1, where we take 9 > 0 so small that R < 1/9 (see Section 2.2 in this 
connection). □ 
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In the course of the proof of Theorem 3.1, which now follows, we usually write 
T — T(t). As mentioned in the Introduction, 


tr(t, D) — \D\p t (0) = / po(t,x,x)dx— / p(t,x,x)dx = — / ro(t,x,x)dx. 


' D 


id 


ID 


We only need to show that 


r D (t, x, x)dx - \dD\C m (t) 


id 


< 


cT 2 
T d R 2 ' 


( 21 ) 


We first consider T = T(t) > R/ 2, and we have 


r D (t,x,x) < 


id 



Pt(0)dx < \D\p t (0) < 4\D\pt(Q) 


T 2 
R 2 ' 


By Lemma 2.5, 


/ °° r \D\ r\D\T 2 ~ d 

r m (t, (q, 0, • • • , 0), (q, 0, • • • , 0 ))dq < -U-T'~ d < 1 l R2 . 

By [5, (23)], we see that (20) holds trivially in this case. 

From now on we assume that T < R/2. For r > 0 we let D r = {x G D : 5 d (x) > r}. 
We have D = D R / 2 U (D\ D R / 2 ). In analyzing the decomposition we shall often use 
our assumptions R < 1/6 and \x — y\ < 1/9, and the heat kernel estimates from 
Lemma 2.6. By Lemma 2.4, 

f DR ro( t ^)d X < C\D Rn \^L_ < CID^ < C|0|^4_. (22) 


Thus, the integral gives insignificant contribution to the trace. 

To handle the integration near dD , we shall estimate the heat remainder of D using 
the heat remainder of halfspace. Let x* G dD be such that \x — x*| = 5d(x). Let 
/ and O be the (inner and outer) balls with radii R such that dl fl 80 = {x*} and 
/ C D C O c . Let lH(a;) denote the halfspace satisfying / C M(x) C O c . By domain 
monotonicity of the heat remainder, and by Lemma 2.2, 


\r D (t,x,x) - r U ( x )(t,x,x)\ < n(t,x,x) - r 0 c(t,x,x) 

= po c (t,x,x ) — p/(t,a:,a:) 

= E x [n < t,x(n ) e O c ;p 0 c(t - t/,x(t/),t)] . 

The next result is an analogue of [1, Proposition 3.1]. 

Proposition 3.2. If T < R/2, then 

E * [r, < UM e O c ; Po <( t - r Mx)] < 1 A . 


The proof of Proposition 3.2 is given in Section 4. 
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Lemma 3.3. IfT< R/2, then 


r D (t,x,x) -r M ( x )(t,x,x) dx 


< 


c\D\T 2 

R 2 T d 


( 23 ) 


' D\D R /2 

Proof. This is an analog of [1, Claim 2] and is proved as follows. By the coarea 
formula and Proposition 3.2 we find that the left side of (23) is bounded above by 


cT r R/2 


RT d 


I dD q 


T d ~^V(T) 

q d - l V(q) 


A 1 ) dq. 


Therefore [1, Corollary 2.14(i)] gives a simplified bound 

c\dD\ f R/2 (T d ~ l V(T) 


A 1 ) dq. 


RT d ~ l J o V q d ~ l V(q) 

The integral over (0, T) is clearly bounded by T. To estimate the integral from T to 
R/2 we note that scaling (9) for q G [T,R/ 2) yields V(T) < C (T/ q)-/ 2 V(q). Also, 

r R/2 


1 -d- 


q 


- /2 dq < / q l - d ~^ 2 dq < oo, 


since d + a/2 > 2. □ 

Recall that r(t, q) = r m (t , (q, 0, • • ■ , 0), (q, 0, • • • , 0)), and C m (t) = f 0 °° r(t, q)dq. 
Lemma 3.4. IfT< R/2, then 


- R/2 


D\D r / 2 

Proof. Using the coarea formula we get 


rw(x)(t,x,x)dx - \dD\ r(t,q)dq 


< 


c\D\T 2 
R 2 T d 


(24) 


D\D R /2 


ru(x)(t, x, x)dx = 


rR/2 


\dD q \r(t,q)dq. 


Hence the left side of the inequality (24) is bounded by 


-R/2 


I \dD q \ — \dD\\r(t, q)dq < 


C\D\ 

R 2 


rR/2 


q r(t, q)dq, 


as follows from [1, Corollary 2.14(iii)]. For q G (0, T\ we have r(t,q ) < pt( 0), hence 


qr(t, q)dq < c 




dq = cT 


2—d 


For the remaining integration, using (17) and (9), we get 


~i/0 


qr(t, q)dq < c 


< c 


-i/0 


t 


< T q d ~ l V 2 {q) 

' 1,9 dq 


dq < c 


- 1/0 


V 2 (T) 
l T q d ~ 1 V 2 {q) 


dq 


-.d -1 — 


< cT 


2—d 




qj q“ 

The last integral converges since d > 2 and a > 0. 


Rlq. 


□ 
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Figure 1. Balls / C D (left), O C D c (right) and P (middle), and 
“a short jump” to point X(ti). Here x E P and |x| = 5j(x). 


Thus, for T < R/2 we have by Lemma 2.4 

c\D\ f 


\dD\ / r(t,q)dq< 
Jr/2 


v2 ^ dq < c ' D 


dq 


CT 2 


R J R /2q d V\q) - R J R/ 2 T d ~ 2 q 2 R?T d ’ 


which is a lower order term. By Lemma 3.3, Lemma 3.4 and (22) we obtain (21). 


4. Proof of Proposition 3.2. 


Let x* = 0, a = (-R, 0,..., 0), b = (R, 0,..., 0), / = B(a,R) and O = B(b,R). 
This also means that x = (t 0 , 0,..., 0) with 0 < x 0 < R/2 , and 5j(x) = |x|, see 
Figure 1. Recall that t < V 2 (R/ 2) or equivalently T < R/2. Before we proceed to 
the heart of the matter we need the following lemma based on spherical integration 
developed in [9, Pages 355-355] and later used in [1, 2], 

Lemma 4.1. For s < R we have 


* (O c \I)CB(0,s) 


dz V(Soc{z)) 
x-z\P P(<M-)) 


< c 


I x\ d+1 ~P/R 

S d +1-P/R 


if /3 > d + 1, 
if (3 < d + 1. 


(25) 


Proof. First we consider V(x) = x a l 2 with a G [0, 2). Let z E A = ( O c \ I) D R(0, s). 
Note that \x — z\ > |x|. If \x — z\ < 2|x|, then \z\ < \x — z\ + |x| < 3|x|, which leads 
to the integral 


' An{|a:—2|<2|a:|} 


dz bo c ( z ) 
x — z\P Sf(z) 


1 f 6&(z) 

X A 8\P J An{\z\<3(\x\As)} $?{z) 
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The last integral is similar to [1, (3.21)]. Using [1, (3.23) and (3.24)] we get the 
following upper bound 

c /-3(NA s) r d _ c (| x | As) d+1 -^ 

\x A s\P J 0 R R 

If \x — z\ > 2|x|, then \x — z\ > \z\/2 and \z\ > \z — x\ — |x| > |x|. By [1, (3.24)], 

f dz S JX1 < c [ 1 s lM dz < £ f r i-K ir , 

>/An{|x—*|>2|*|} 1*^ ^1^ ^1 (^0 ./An{s>p|>|a;|} ^7 (") ~^./|a;|As 

If /3 > d + 1, then the last integral is bounded by c|a;| d+1_/3 , while for (5 < d + 1 we 
get the upper bound cs^ 1 ^ 13 . 

This settles (25) for V(x) = x a ^ 2 with a £ [0,2). Note that the form of the right 
hand side of (25) does not depend on a. 

Consider general -0 £WUSC(a) and the corresponding ladder-height function V. 
Due to the scaling property (10) we have 


V(Sqc(z)) hg c (z) 

V{5 t (z)) ~ ' 8f(z) ’ 


if 8o c (z) > 8i{z). 


If 5o c (z) < Si(z), then the fraction is bounded by 1, since V is monotone. Therefore, 
we can use the previous special case with a = a and a = 0 to finish the proof. □ 


We return to the core proof of Proposition 3.2. In view of Lemma 2.3 we want to 
estimate 


E' r [t/ < t, X (tj) £ O c \p 0 c(t - n, X (r/), x)] 


Pi(s,x,y) / v(y — z)p 0 c(t 


'I Jo 

= h + h + h-, 


'o c \i 


s, x , z) dzdsdy 


which splits the integration into three subregions, as specified and estimated below: 
h : \z\ > Rj 2, 

/ 2 : t/2 < s < t AND \x — z\ < T AND \z\ < Rj 2, 

R : (s < t/2 OR \x — z\ > T) AND \z\ < Rj2. 

The setting, especially that of J 2 , is illustrated on Figure 1. 


4.1. Long jump: integral I\. On I\ we have \z\ > R/2, hence \x — z\ > R/ 3, thus 
by (13) 


h = 


< 


1 Jo 


Pi{s,x,y) 

ct 


R d V 2 (R/ 3) 
ct 

R d V 2 (R/ 3) 


I v(y — z)p{t — s, z, x)dsdzdy 

\z\>R/2 

pi(s,x,y ) / u(y — z)dsdzdy 
J p c 

cV 2 {T) 


¥ x (r I <t,\X(r I )\>R/2)< 


R d V 2 (R/2y 
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where the last inequality follows from sublinearity (8) of V. Since T < R/2, we have 

cV 2 {T) c c 

R d V 2 {R/2) -R*~ RT d - l ' 

Since |x| < R/2, by monotonicity of V we get 

cV 2 (T) cV(T) cV(T) 

R d V 2 (R/2) - R d V(R/2) ~ fi|z|'‘->r(|z|)' 


4.2. Long exit time and short jump: integral / 2 . Here we have |x| < \x — z\ < T, 

and \z\ < \x — z\ + |x| < 2 T. By Lemma 2.6, t/2 < q < T and (12), 


Pi(q,x,y) < T 


V(T) ' 


Let S = ( O c \/)n{|z| <2T}. We get the following upper bound, 


h = 


P!(q,x,y) / u(y - z)p 0 c(t - q,z,x)dqdzdy 


I Jt/2 


< c J T(t) 
cT~ d 




< 



V(T) J s \y-z\*V*(\y-z\) 
V(Sj(z)) G' 0 l (rtm (x,z 


V(T) J S J, \v - z\ d V(\y - z\) V(S,(z)) 

where we use Sj(y) < \y — z\. Scaling (9) gives 


(x, z)dzdy 

dydz, 


h < 


cT 


Sf /2 (z) Gt {m \x,z) 
V(T) J s \y - z\ d +zl 2 V(S,(z)) 



dydz. 


We then rewrite the inner integral in spherical coordinates, use Green function esti¬ 
mate (18) and |x| < T , 


I2 < 


cT 


—d roo X—R 


_ ef‘{z)dr [ V(\x\)V(6 0 '(z)) 

VC T)J sm ri+s/2 j s \x-z\*V(6,(z)) 


dz 


< cT 


-d 


dr 


V(S 0 c(z)) 


dz = cT 


f x r 1+ -/ 2 Js \x - z\ d V(5 I (z)) 
Using Lemma 4.1 with f3 = d and s = 2T we get 

cT i-d 


- d [ V(Sqc(z)) 

l s \x - z^Vid^z)) 


dz. (26) 


I2 < 


R 


Since |x| < T, we get the desired estimate from Proposition 3.2. 
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4.3. Short exit time or medium jump: integral J 3 . Let S = ( 0 C \ /) D {|z| < 
R/ 2}. We have \x — z\ > T or s < t/2. In either case, Lemma 2.6 and sublinearity of 
V implies 


Po c (t — s,x, z) < [T d A 


V 2 (T) 


V{8 0c {z)) 


Therefore by Lemma 2.9, 

- V\R/2) r- 

Pi(s,x,y ) / u(y - z) T~ d A 


h< 


= c 


< c 


I J 0 Js 

J kJ 2(R) {x, z ) (T~ d A 


x - z\ d V 2 (\x - z\) J V(T) ' 

V\T) \ V(SoM) 


x - z\ d V 2 (\x - z\) J V(T) 
V 2 (T) ^ V(5 0e (^)) 


dzdsdy 


-dz 


v(M) 


<S V(5!(z)) \x - z\ d 


T~ d A 


z\ d V 2 {\x - z\) J V(T) 

V 2 (T) \ V(S 0e {z)) 

V(T) 


dz. 


\x — z\ d V 2 (\x — z\) 

If \x — z\ < T, then we are satisfied with T~ d from the minimum and we note 
IA(|:r|) < V(T). We arrive at (26), and hnish the proof in the same way as in the 
previous cases. 

We are left with the case \x — z\ > T, and we have 

V(M) V(5 0c (z)) 


h < cV{T ) 

Since ifj gWLSC(o), we get 
h<cV{T) r 


l s \x — z\ 2d V 2 (\x — z\) V(5i(z)) 

|s/2 nsoM) 


-dz. 


\x\ 


s \x — z\ 2d+ -/ 2 V[\x — z\) V(5i(z)) 


dz 


< 


<W(T)|:r|-/ 2 


v(SoM) 


-dz, 


(:T V \x\) d ~ l V{T V |x|) Js \x - z\ d+1 +^ 2 v(6j(z)) 

where the last inequality follows from the monotonicity of V, since \x — z\ > |x| V T. 
Now we use Lemma 4.1 with (3 = d + 1 + a/2, to get 

/,<_4W_, 

- (T V |x|) d_1 r(T V |x|)i? 

Here the right hand side is comparable with the required upper bound. 
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